							
							
							
*---- Authors: Maria Nordbrandt, Gina Gustavsson, Karen Nielsen Breidahl
*---- The Unifying Magic, EJPR 2025


*This script contains code used to produce the main results reported in Figure 1 in the main text.
*Code to produce Figure 1 can be found in the R script called Figure 1 and D1.R.  



							**FIGURE 1
										
	
									*Sweden
									*Flag
										eststo clear
										eststo: reg marry_n treat_flag  if id_se!=. 
										
										
										eststo est1 
										eststo: reg marry_n_w treat_flag  if id_se!=.
										
										
										eststo est2 
										eststo: reg marry_dist_n treat_flag  if id_se!=.  
										
										
										eststo est3
										eststo: reg marry_dist_n_w treat_flag  if id_se!=.
										eststo est4 
										
										eststo: reg dislike_n treat_flag  if id_se!=.  
										
										
										eststo est5
										eststo: reg dislike_n_w treat_flag  if id_se!=.
										
										
										eststo est6
										eststo: reg dislike_dist_n treat_flag  if id_se!=.  
										
										
										eststo est7
										eststo: reg dislike_dist_n_w treat_flag  if id_se!=.
										eststo est8
										
										eststo: reg traits_n treat_flag  if id_se!=.  
										
										
										eststo est9
										eststo: reg traits_n_w treat_flag  if id_se!=.
										
										
										eststo est10
										eststo: reg traits_dist_n treat_flag  if id_se!=.  
										
										
										eststo est11
										eststo: reg traits_dist_n_w treat_flag  if id_se!=.
										eststo est12
										
										esttab, replace ///
se stats(r2_a  N, fmt(a2) labels("Adj. R2" "Observations")) ///
b(%9.3fc)  compress star(* 0.10 ** 0.05 *** 0.01) /// 
unstack label nobaselevels nomtitles eqlabels(none) obslast nonotes  //////
varlabels(_cons Constant)
										
										*Cake
										eststo clear
										eststo: reg marry_n treat_cake  if id_se!=.  
										
										
										eststo est1 
										eststo: reg marry_n_w treat_cake  if id_se!=.
										
										
										eststo est2 
										eststo: reg marry_dist_n treat_cake  if id_se!=.  
										
										
										eststo est3
										eststo: reg marry_dist_n_w treat_cake  if id_se!=.
										eststo est4 
										
										eststo: reg dislike_n treat_cake  if id_se!=.  
										
										
										eststo est5
										eststo: reg dislike_n_w treat_cake  if id_se!=.
										
										
										eststo est6
										eststo: reg dislike_dist_n treat_cake  if id_se!=.  
										
										
										eststo est7
										eststo: reg dislike_dist_n_w treat_cake  if id_se!=.
										eststo est8
										
										eststo: reg traits_n treat_cake  if id_se!=.  
										
										
										eststo est9
										eststo: reg traits_n_w treat_cake  if id_se!=.
										
										
										eststo est10
										eststo: reg traits_dist_n treat_cake  if id_se!=.  
										
										
										eststo est11
										eststo: reg traits_dist_n_w treat_cake  if id_se!=.
										eststo est12
										
										esttab, replace  ///
se stats(r2_a  N, fmt(a2) labels("Adj. R2" "Observations")) ///
b(%9.3fc)  compress star(* 0.10 ** 0.05 *** 0.01) /// 
unstack label nobaselevels nomtitles eqlabels(none) obslast nonotes  //////
varlabels(_cons Constant)
										
										*Denmark
										*flag
										eststo clear
										eststo: reg marry_n treat_flag  if id_dk!=.  
										
										
										eststo est1 
										eststo: reg marry_n_w treat_flag  if id_dk!=.
										
										
										eststo est2 
										eststo: reg marry_dist_n treat_flag  if id_dk!=.  
										
										
										eststo est3
										eststo: reg marry_dist_n_w treat_flag  if id_dk!=.
										eststo est4 
										
										eststo: reg dislike_n treat_flag  if id_dk!=.  
										
										
										eststo est5
										eststo: reg dislike_n_w treat_flag  if id_dk!=.
										
										
										eststo est6
										eststo: reg dislike_dist_n treat_flag  if id_dk!=.  
										
										
										eststo est7
										eststo: reg dislike_dist_n_w treat_flag  if id_dk!=.
										eststo est8
										
										eststo: reg traits_n treat_flag  if id_dk!=.  
										
										
										eststo est9
										eststo: reg traits_n_w treat_flag  if id_dk!=.
										
										
										eststo est10
										eststo: reg traits_dist_n treat_flag  if id_dk!=.  
										
										
										eststo est11
										eststo: reg traits_dist_n_w treat_flag  if id_dk!=.
										eststo est12
										
										esttab, replace  ///
se stats(r2_a  N, fmt(a2) labels("Adj. R2" "Observations")) ///
b(%9.3fc)  compress star(* 0.10 ** 0.05 *** 0.01) /// 
unstack label nobaselevels nomtitles eqlabels(none) obslast nonotes  //////
varlabels(_cons Constant)
									
										
										*cake
										eststo clear
										eststo: reg marry_n treat_cake  if id_dk!=.  
										
										
										eststo est1 
										eststo: reg marry_n_w treat_cake  if id_dk!=.
										
										
										eststo est2 
										eststo: reg marry_dist_n treat_cake  if id_dk!=.  
										
										
										eststo est3
										eststo: reg marry_dist_n_w treat_cake  if id_dk!=.
										eststo est4 
										
										eststo: reg dislike_n treat_cake  if id_dk!=.  
										
										
										eststo est5
										eststo: reg dislike_n_w treat_cake  if id_dk!=.
										
										
										eststo est6
										eststo: reg dislike_dist_n treat_cake  if id_dk!=.  
										
										
										eststo est7
										eststo: reg dislike_dist_n_w treat_cake  if id_dk!=.
										eststo est8
										
										eststo: reg traits_n treat_cake  if id_dk!=.  
										
										
										eststo est9
										eststo: reg traits_n_w treat_cake  if id_dk!=.
										
										
										eststo est10
										eststo: reg traits_dist_n treat_cake  if id_dk!=.  
										
										
										eststo est11
										eststo: reg traits_dist_n_w treat_cake  if id_dk!=.
										eststo est12
										
										esttab, replace ///
se stats(r2_a  N, fmt(a2) labels("Adj. R2" "Observations")) ///
b(%9.3fc)  compress star(* 0.10 ** 0.05 *** 0.01) /// 
unstack label nobaselevels nomtitles eqlabels(none) obslast nonotes  //////
varlabels(_cons Constant)

